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Abstract 

Time dependent quadratic Hamiltonians are well known as well in classi- 
cal mechanics and in quantum mechanics. In particular for them the corre- 
spondance between classical and quantum mechanics is exact. But explicit 
formulas are non trivial (like the Mehler formula). Moreover, a good knowlege 
of quadratic Hamiltonians is very useful in the study of more general quan- 
tum Hamiltonians and associated Schrodinger equations in the semiclassical 
regime. 

Our goal here is to give our own presentation of this important subject. We 
put emphasis on computations with Gaussian coherent states. Our main mo- 
tivation to do that is application concerning revivals and Loschmidt echo. 
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1 Introduction 

This paper is a survey concerning exact useful formulas for time dependent Schrodinger 
equations with quadratic Hamiltonians in the phase space. One of our motivations 
is to give a detailed proof for the computation of the Weyl symbol of the propagator. 
This formula was used recently by Melhig- Wilkinson to suggest a simpler proof 
of the Gutzwiller trace formula [Sj. 

There exist many papers concerning quantum quadratic Hamiltonians and exact 
formulas. In 1926, Schrodinger [21] has already remarked that quantification of the 
harmonic (or Planck) oscillator is exact. 

The best known result in this field is certainly the Melher formula for the harmonic 
oscillator (see for example [5]). 

Quadratic Hamiltonians are very important in partial differential equations on 
one side because they give non trivial examples of wave propagation phenomena and 
in quantum mechanics and on the otherside the propagation of coherent states by 
general classes of Hamiltonians, including —h 2 A + V, can be approximate modulo 
0{h°°) by evolutions of quadratic time dependent Hamiltonians [Tj ITT^ ITij. 

In his works on pseudodifferential calculus, A. Unterberger in [221 EHj has given 
several explicit formulas connecting harmonic oscillators, Gaussian functions and 
the symplectic group. This subject was also studied in [HI HHJ E2| • More recently de 
Gosson P has given a different approch for a rigorous proof of the Melhig- Wilkinson 
formula for metaplectic operators, using his previous works on symplectic geometry 
and the metaplectic group. 

Here we shall emphasis on time dependent quadratic Schrodinger equation and Gaus- 
sian Coherent States. It is well known and clear that this approach is in the heart 
of the subject and was more or less present in all papers on quantum quadratic 
Hamiltonians. In this survey, we want to give our own presentation of the subject 
and cover most of results appearing in particular [T31 E] . 

Our main motivation to revisit this subject was to prepare useful tools for applica- 
tions to revivals and quantum Loschmidt echo [2] . We shall see that in our approach 
computations are rather natural direct and explicit. 

2 Weyl quantization. Facts and Notations 

Let us first recall some well known facts concerning Weyl quantization (for more 
details see HUM. 

The Planck constant h > is fixed (it is enough to assume K = 1 in the homogeneous 
quadratic case). 

The Weyl quantization is a continuous linear map, denoted by • or by Op w , de- 
fined on the temperate Schwartz space distribution iS'(IR 2n ) into £ w (S(R n ), S'(R n )) 
where C W (E,F) denotes the linear space of continuous linear map from the linear 
topological E into the linear topological F with the weak topology. 
M 2 ™ is a symplectic linear space with the canonical symplectic form a(X, Y) = JX-Y 
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where J is defined by 




and t n denotes the identity n x n matrix. 

Let us introduce the symplectic group Sp(2n): it is the set of linear transformations 

of R 2n preserving the 2-form a. 

For X G M 2n , we denote X = (q,p) eK"x R n . 

The Weyl quantization is uniquely determined by the following conditions: 
(WO) A i— > A is continuous. 

(Wl) (#)(g) = #(g), (W)(q) = D q ^{q), (D q = ^) 
for every ip E S(W n ). 

(W2) exp[i(a -q+~P-p)] = exp[i(a ■ q + (3 ■ p)], V a,(3eW l 

Let us remark that a ■ q + (3 ■ p is self-adjoint on L 2 (IR n ) so, exp[i(a ■ q + (3 ■ p)] 
is unitary. In particular, if z = (x,£) then 

T(z) =: exp[-(£-g-x-p) 

is the quantized translation by z in the phase space (Weyl operators). 
From (Wl), (W2), using continuity and H- Fourier transform, defined by 

A(Y) = ( e-r> x - Y A{X)dX 

we have 

A=(2nh)- 2n [ A(a, (3) exp{Ua ■ q + (3 ■ p)}dadp (2.1) 

In general, equality (|2.1|) is only defined in a weak sense i.e through the duality 
bracket between S' and S, < A(p,ip > for arbitrary (p, ip e tS(R n ). 

Definition 2.1 A zs i/ie Weyl contravariant symbol of A if they are related through 
the formula \2.1)) . 

Using the explicit action of Weyl operators on iS(M n ) and Fourier analysis, we 
get the following formula (see JUICE]) 

(A<p)(x) = (2nh)- n Jl^expl^x-y) • £L4 (^-,t) ^)«- (2-2) 
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In particular, if Ka denotes the Schwartz kernel of A, we have 

f % { u u\ 

M<l,v) = J exp[--u ■ p]K A [q+ -,q- -J du (2.3) 

These formulas are true in the distribution sense in general, and pointwise if A is 
smoothing enough ( surely if for example A G iS(K n )). 

Let us remark that they are consistent and that A i— > A is a bijection from iS'(M 2n ) 
into £ w (S(M? n ), iS'(R 2n )). In particular we have the following inversion formula 

Proposition 2.2 For every A G C w (S(M. 2n ),S'(M. 2n )), there exists a unique con- 
travariant Weyl symbol A G 5'(R 2n ) given by the following formula 

A{X) = 2 n Tr[AS ym (X)} (2.4) 

where S ym (X) is the unitary operator in L 2 (IR n ) defined by 

S ym {X)ip{q) = (7r/r n e^ (a; ~ 9 V(2x - q) 

for X = (x, f ) . 

Sketch of proof: 

We first prove the formula for A G C w (S'(R 2n ), S (K 2n )) .The general case follows 
by duality and density. 

We start with the following formula, easy to prove if A, B G 5(M 2n ), 

Tr[AB] = (2nh)- n [ A{X)B{X)dX (2.5) 

Assume now that B is a bounded operator in L 2 (R n ). The Weyl symbol B of B 
satisfies: 

Tr[AB] = (2nh)- n < A, B > 1 (2.6) 

we have to check that the Weyl symbol of S ym (X) is {nh) n 5x where 5x is the Dirac 
mass in X. 

To prove that let us consider any ip G iS(M n ) and denote by W v the Weyl symbol 
of the projector ip i— >< ip, ip > if. W v is called the Wigner function of if. A direct 
computation using (|2.3|) gives 

ft u u 

W v (q,p)= / exp[--u-p](p{q+-)<p{q- -)du 

We get the proposition applying formula (j2.6|) with B = S ym (X) and A = tc^.D 

We shall see later that it may be convenient to introduce the covariant Weyl 
symbol for A which has a nice connection with Weyl translations. 

1 the bracket <,> denotes the usual bilinear form (integral or distribution pairing). We shall 
denote (•!•) the Hermitean scsquilinear form on Hilbert spaces, linear in the second argument. 



5 



Proposition 2.3 For every A e C w (S(M. 2n ), iS'(M 2n )) there exists a unique temper- 
ate distribution A# on M? n , named covariant Weyl symbol of A, such that 

A = (2nh)- n [ A*{X)f{X)dX (2.7) 

Moreover we have the inverse formula 

A*{X) = Tr[Af(-X)} (2.8) 

( As above, if A is not trace-class, this formula has to be interpreted in a weak 
distribution sense). 

The covariant and contravariant Weyl symbols are related with the following formula 

A*(X) = (2nhy n A(JX). (2.9) 
A o J is named the symplectic Fourier transform of A 

Proof 

These properties are not difficult to prove following for example [20J.D 



We define the (usual) Gaussian coherent states tp z as follows: 

(p z = f(z)<po, VzGt 2n (2.10) 

where 

<p (x) := (nh)' n/4 e~ x2/2h (2.11) 
We get the following useful formula for the mean value of observables: 

Corollary 2.4 With the above notations, for every <p, ip e «S(R n ), we have 

(<p\Ai/>) = (2nh)- n [ A*(X)(<p\f(X)i))dX (2.12) 
In particular for Gaussian Coherent States, we have 

( Vz \A Vo ) = (2nK)- n A*(X) exp ( "^^ " £f(X, «)) dX (2.13) 
Proof: 

The first formula is a direct consequence of definition for covariant symbols. 

The second formula is a consequence of the first and the following easy to prove 

equalities 

f(z)f(z') = e W (^(z,z')^f(z + z') (2.14) 

{<Pz\<Po) = e « (2.15) 

□ 

For later use let us recall the following 
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Definition 2.5 Let be A, B G C w (S{R n ), <S'(M n )) such that the operator compo- 
sition AB is well defined. Then the Moyal product of A and B is the unique 
A#B G S'(R 2n ) such that 

AB = A#B (2.16) 
For details concerning computations rules and properties of Moyal products see 

huhbi 



3 Time evolution of Quadratic Hamiltonians 

In this section we consider a quadratic time-dependent Hamiltonian, 
Ht{z) = J2i<j k<2n c j,k(t)zjZk, with real and continuous coefficients Cj t f-(t), defined 
on the whole real line for simplicity. It is convenient to consider the symplectic 
splitting z = (q,p) 6 I" x 1" and to write down H t (z) as 

H t (q,p) = ^ ( G tQ ■ Q + 2L tQ ■ P + K t p ■ p) 

where K t ,L t , G t are real n x n matrices, K t and G t being symmetric. 
The classical motion in the phase space is given by the linear equation 

L J)(l), (3-17) 



p J ' \ L t K t J \p 

where L T is the transposed matrix of L. This equation defines a flow, F t (linear 
symplectic transformations) such that F = t. On the quantum side, H t is a family 
of self-adjoint operators on the Hilbert space 7i = L 2 (M. n ) (this will be proved 
later). The quantum evolution follows the Schrodinger equation, starting with an 
initial state G Ti. 

in ~tt = Ht ^ U = ^ < ^ 3 ' 18 ' ) 
Suppose that we have proved existence and uniqueness for solution of (j3.18|) . we 
write ip t = U t <p. The correspondence between the classical evolution and quantum 
evolution is exact. For every A G 5(IR 2 ' 1 ), we have 

Proposition 3.1 

Ut.A.Uf 1 = AF t (3.19) 

Sketch of proof 

For any quadratic Weyl symbol B we have the exact formula 

^iB,A] = {B~A} (3.20) 

where {B, A} = VB ■ JV A denote the Poisson bracket, and [B, A] = B.A - A.B is 
the Moyal bracket. So we can prove Proposition (|3.1j) by taking derivative in time 
and using ()3.20|) (see [IH] for more details). | 
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Now we want to compute explicitly the quantum propagator U to j in terms of 
classical evolution of H t . 

One approach is to compute the time evolution of Gaussian coherent sta(es, Ut 0t tfz, 
or in other word to solve the Schrodinger equation (j3.18j) with <p = <p z , the Gaus- 
sian coherent state in z G M". Let us recall that (p z = T(z)ipo and <po(x) = 

(7T^)- n / 4 exp 



2fi 



4 Time evolution of Coherent States 

The coherent states system {ip z } z ^R2 n introduced before is a very convenient tool 
to analyze properties of operators in L 2 (R n ) and their Schwartz distribution kernel. 
To understand that let us underline the following consequence of the Plancherel 
Formula for the Fourier transform. In all this section we assume h = 1. For every 
u E L 2 (R n ) we have 

\u(x)\ 2 dx = (2ir)- n [ \(u,<p z )\ 2 dz (4.21) 
Jul 2 " 

Let R be some continuous linear opertor from «S(M n ) into S'(M. n ) and Kr its Schwartz 
distribution kernel. By an easy computation, we get the following representation 
formula 

K R (x,y) = (2tt)"" f (R Vz )(x)^Mdz. (4.22) 

In other words we have the following continuous resolution of the Schwartz distri- 
bution kernel of the identity 



5(x - y) = (2tt) n / ip z {x)ip z (y)dz. 

JR 2n 

This formula explains why the Gaussian coherent system may be an efficient tool 
for analysis of operators on the Euclidean space M n . 
Let us consider first the harmonic oscillator 

It is well known that for t ^ kir, fceZ the quantum propagator e~ ttH has an explicit 
Schwartz kernel K(t;x,y) (Mehler formula). 

It is easier to compute with the coherent states ip z . (p is an eigenstate of H, so we 
have 

e-^o = e- l4/ Vo (4.24) 
Let us compute e~ ltH ip z , Vz G M 2 , with the following ansatz 

e-* 6 ^ = e iSt ^f (2t)e _it/ Vo (4-25) 
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where z t = (qt,Pt) is the generic point on the classical trajectory (a circle here), 
coming from z at time t = 0. Let be ipt,z the state equal to the r.h.s in (|4.25jl . and 
let us compute S t (z) such that ip t ,z satisfies the equation i^ip = Hip <p\t=o = i>o,z- 
We have 

f{z t )u{x) = e l{ptX - qm/2) u{x - q t ) 

and 

= e i{Hz)-t/2+ Pt x- qm /2)^ x _ ^ ( A 2 ^ 

So, after some computations left to the reader, using properties of the classical 
trajectories 

qt =Pt Pt = ~qt, Pt + qt = P 2 + Q 2 , 

the equation 

^t, z (x) = \( D l + x 2 )*PtAx) (4-27) 

is satisfied if and only if 

bt{z) = \{p t qt-pq) (4.28) 

□ 

Let us now introduce the following general notations for later use. 
F t is the classical flow with initial time to = and final time t. It is represented as 
a 2n x In matrix which can be written as four n x n blocks : 

C t D, ) ■ < 4 ' 29 > 
Let us introduce the following squeezed states. ip r is defined as follows. 

ip T {x) = a r exp ( -^—Tx ■ x ) (4.30) 
\2ft / 

where r G S n , S„ is the Siegel space of complex, symmetric matrices T such that 
S(r) is positive and non degenerate and ar G C is such that the L 2 -norm of ip r is 
one. 

We also denote ip\ = T(z)ip r . 
For r = it, we denote ip = ip l1 . 

Theorem 4.1 We have the following formulae, for every x G M™ and z G M? n , 

U t <p r (x) = <p Tt (x) (4.31) 
U t <fi(x) = f(F t z)<p Vt (x) (4.32) 

where T t = (C t + iD t T)(A t + iB t T)- 1 and a Tf = a r (det(A t + iB t T)y 1/2 . 
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Beginning of the proof 

The first formula can be proven by the ansatz 

U t (p (x) = a(t) exp (j^ t x ■ x 

where r t e S n and a(t) is a complex values time dependent function. We get first 
a Riccati equation to compute r t and a linear equation to compute a(t). 
The second formula is easy to prove from the first, using the Weyl translation oper- 
ators and the following known property 

U t T(z)U; = T(F t z). 

Let us now give the details of the proof for z = 0. 

We begin by computing the action of a quadratic Hamiltonian on a Gaussian (h = 
!)• 



Lemma 4.2 



Lx ■ D x ei Tx - x = (L T x ■ Tx - -TrL)e5 r " 

2 



Proof 

This is a straightforward computation, using 



Lx ' D x — — ^ ^ Lj 

and, for uj G W 1 , 



XjD k + D k xj 

l<j,k<n 



{uj.D x )e^ x - x = (Tx ■ u)e^ x - x 

□ 

Lemma 4.3 

(GD X ■ D x )^ Vxx = (GTx ■ Tx - iTr(GT)) ^ Tx ' x 

Proof 

As above, we get 

H^ Txx = (]-Kx -x + x-LTx + ]-GYx ■ Tx - ^Tr(L + GT)] e^™ (4.33) 
We are now ready to solve the equation 

= Hi, (4.34) 

with 

M^oix) = g(x) := (2n)-^e- x2 / 2 . 
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We try the ansatz 

#,i) = a(t)e^ r '" (4.35) 

which gives the equations 

t t = -K-2T?L-T t GT t (4.36) 

/(*) = -\(Tr{L + GT t ))f(t) (4.37) 

with the initial conditions 

T = it, a(0) = (2vr)- n / 2 

Remark: T T L et LT determine the same quadratic forms. So the first equation is 
a Ricatti equation and can be written as 

f t = -K- Y t L T - LT t - T t GT t , (4.38) 

where L T denotes the transposed matrix for L. We shall now see that equation 
()4.38|) can be solved using Hamilton equation 

F t = J ( fr n) F t (4-39) 



F = 1 (4.40) 



We know that 

' A t B t 



C t D t 

is a symplectic matrix Vt. So we have det(A t + iB t ) ^ Wt (see below). Let us 
denote 

M t = A t + iB t , N t = C t + iD t (4.41) 
We shall prove that T t = NtMf 1 . By an easy computation, we get 



M t = L T M t + GNt 



N t = -KMt-LNt (4.42) 



Now, compute 
d 



(NtMr 1 ) = NM~ l - XM '.\7.\/ 1 
at 

= -K - L.XU 1 - \ 1 (l.'M + G.\).\l : 

= -K - I..XM 1 - XM '/.' - NM~ 1 GNM' 1 (4.43) 

which is exactly equation ()4.38|) . 

Now we compute a(t), using the following equality, 

Tr (L T + G(C + iD){A + iB)' 1 ) = Tr(M)M _1 = Tr (L + GT t ) 
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using TrL = TrL T . Let us recall the Liouville formula 

^ log(det M t ) = Tr^M- 1 ) (4.44) 
at 

which give directly 

a(t) = (27t)-™/ 2 (det(A + iB t )Y 1/2 (4.45) 
To complete the proof, we need to prove the following 

Lemma 4.4 Let S be a symplectic matrix. 

A B 
C D 

Then det(A + iB) ^ and 3(C + iD)(A + iB)- 1 is positive definite. 

We shall prove a more general result concerning the Siegel space E„. 
Lemma 4.5 If 



S = 



A B 
C D 



is a symplectic matrix and Z G E„ then A+BZ et C+DZ are non singular and 
(C + DZ){A + BZ)- 1 E E n 

Proof 

Let us denote E := A+BZ, F := C+DZ. F is symplectic, so we have F T JF = J. 
Using 

we get 

(f,F T )j(^(/,Z)j(^0 (4.46) 

which gives 

E T F = F T E 

In the same way, we have 

f) = I(/, (4.47) 

We get the following equation 

F T £ - £ T F = 2i%Z (4.48) 
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Because is non degenerate, from (|5.68jl . we get that E and F are injective. If 

x G C™, Ex = 0, we have 

Ex = x T E T = 

hence 

x T QZx = 

then x = 0. 

So, we can define, 

a(S)Z= (C + DZ)(A + BZ)- 1 (4.49) 
Let us prove that a(S) G £„. We have: 

a(M)Z = (a(M)Z) T = (E^f F T = (E^ 1 ) 7 E T FE' 1 = EE' 1 = a(M)Z. 

We have also: 

tFE' 1 — EE' 1 - F T E — E T F _ 

E T E = = QZ 

2% 2i 

and this proves that &(a(M)) is positive and non degenerate. 
This finishes the proof of the Theorem for z — . □ 

Remark 4.6 For a different proof of formula \4-31 ), using the usual approach of 
the metaplectic group, see the book J<^. 

The family {(p z } z m 2n spans all of L 2 {R n ) (see for exemple JTffi for properties of 
the Fourier- Bargmann transform) so formula \4.31\) wholly determines the unitary 
group Ut- In particular it results that Ut is a unitary operator and that H t has a 
unique self-adjoint extension in L 2 (M. n ). This is left as exercises for the reader. 

Remark 4.7 The map S i— > a(S) defines a representation of the symplectic group 
Sp(2n) in the Siegel space E n . It is easy to prove that a(SiS 2 ) = a(Si)a(S 2 ) ■ This 
representation is transitive. Many other properties of this representation are studied 
in JEy. 



5 The metaplectic group and Weyl symbols com- 
putation 

A metaplectic transformation associated with a linear symplectic tranformation F G 
Sp(2n) in R 2n , is a unitary operator R(F) in L 2 (R n ) satisfying one of the following 
equivalent conditions 

R(F)*AR(F) = VA e 5(l 2n ) (5.50) 

R(F)*T(X)R(F) = T[F~ 1 (X)], VX G M 2n (5.51) 
R(F)*AR(F) = 5UF, 
for A(q,p) = qj, 1 < j < n and A(q,p) — Pk, 1 < k < n. (5.52) 
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We shall see below that for every F £ Sp(2n) there exists a metaplectic transforma- 
tion R(F). 

Let us remark that if R\(F) and R,2(F) are two metaplectic operators associ- 
ated to the same symplectic map F then there exists A £ C, |A| = 1, such that 
Ri(F) = XR 2 (F). It is also required that F \— > R(F) defines a projective represen- 
tation of the real symplectic group Sp(2n) with sign indetermination only. More 
precisely, let us denote by Mp(n) the group of metaplectic transformations and 
7T P the natural projection: M p — > Sp(2n) then the metaplectic representation is 
a group homomorphism F \— > R(F), from Sp(2n) onto Mp(n)/{1, — 1}, such that 
7Tp[R(F)] = F, WF £ Sp(2n) (for more details concerning the metaplectic represen- 
tation see We shall show here that this can be achieved straightforward using 
Theorem 14.11 

For every F £ Sp(2n) we can find a C 1 - smooth curve F t , t £ [0, 1], in Sp(2n), 
such that Fq = t and F\ = F. An explicit way to do that is to use the polar 
decomposition of F, F = V\F\ where V is a symplectic orthogonal matrix and 
\F\ = a/ F t F is positive symplectic matrix. Each of these matrices have a logarithm, 
so F = e K e L with K, L Hamiltonian matrices, and we can choose F t = e tK e tL . Any 
way, Ft is clearly the linear flow defined by the quadratic Hamiltonian H t (z) = ~SfZ-z 

where S t = —JFfFf 1 . So using above results, we define R(F) = U\. From 
this definition and Theorem 14.11 we can easily recover the usual properties of the 
metaplectic representation. 

Proposition 5.1 Let us consider two symplectic paths F t and F[ joining 1 (t = 0) 
to F (t — 1). Then we have U\ = ±U{ (with obvious notations). 
Moreover, if F 1 ^ 2 £ Sp(2n) then we have 

R{F l )R(F 2 ) = ±R(F 1 F 2 ). 

Proof 

Using (J4.31j) we see that the phase shift between the two paths comes from variation 
of argument between and 1 of the complex numbers b(t) = det(A t + iB t ) and 
b'(t) = det(A' t + iB' t ). 

We have arg[fo(£)] = S ^ J ' jj^ds J and it is well known (see Lemma (|5.4jl below and 
its proof) that 

with N £ Z. So we get 

6(1)-V2 = e ^, (1) -l/2 

The second part of the proposition is an easy consequence of Theorem 14.11 concern- 
ing propagation of squeezed coherent states with little computations. □ 
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In a recent paper ^\ the authors use a nice explicit formula for the Weyl symbol 
of metaplectic operators R{F). In what follows we detail a rigorous proof of this 
formula including computation of the phase factor. In principle we could use Theo- 
rem H]T] to compute the Weyl symbol of the propagator U t0it - But in this approach 
it seems difficult to compute phase factors (Maslov- Conley-Zehnder index). 
For technical reason, It is easier for us to compute first the contravariant Weyl sym- 
bol, U t , for the propagator U t defined by H t . In any case, U t is a Schwartz temperate 
distribution on the phase space R 2n . 
We follow the approach used in Fedosov [7j. 

It is enough to assume h = 1 . In a first step we shall solve the following problem 

^|p = m 

Ul = t e (5.53) 

where t e is a smoothing family of operators such that lim t e = 1. It will be conve- 
nient to take l e (X) = exp(— £:|X| 2 ). 



Let us recall that # denotes the Moyal product for Weyl symbols. So for the 
contravariant symbol U t (X) of U t we have 

i^U t (X) = (H t #U t )(X) (5-54) 
Because H t is a quadratic polynomial we have 

{H t 4U t ){X) = H t (X)U t (X) + l{# t , U t }(X) 

1% 

-\(p& v - d x d v ) 2 H t (X)U t (Y)\ x=Y (5.55) 

o 

where d x = X = (x,£), Y = (y,r}). 

It seems natural to make the following ansatz 

U t (X) = a(t)E t (X), where 

E t (X) = exp(iM t X-X). (5.56) 

a(t) is a complex time dependent function, M t is a time dependent 2n x 2n complex, 
symmetric matrix such that ^sM t is positive and non degenerate. 
A, B being two classical observables, we have: 

{A , B} = VA . JVB 

and 

(d x d v -d y ds) 2 A(x,0B(y,rj)\ x=Y = d 2 x2 A d^B + d 2 x2 B d 2 e A - 23%A 8%B 

= -Tr(JA"JB") (5.57) 
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where A" is the Hessian of A (and similarly for B). Applying this with 



we get: 



A(X) = ^S t X.X, B(X) = exp(iM t X.X) 



H t #E t (X) = H t (X) E t (X) + JS t X . M t X E t {X) + l -Ti(JS t JB") 



However, 

VB = 2iB(X) M t X 
{B") jh = 2% ((M t ) jk + 2i{M t X) j (M t X) k ) B(X) 

so that: 

^Tr(JS t JB") = (~Tv(JS t JM t ) - l -M t X . JS t JM t Xj B(X) 
Therefore the Ansatz (|5.56|) leads to the equation: 

ia(t) - a{t)M t X . X =~ (S t X . X + M t JS t X . X - S t JM t X . X) a{t) 

+ t -a(t)Tr(MtS t ) - l -a(t)M t X . S t M t X (5.58) 

where we have introduced the Hamiltonian matrices 

M t := JM t , S t := JS t 
Then equation ()5.58j) is equivalent to 

M t = ^(M t + l)S t (M t -l) (5.59) 
a t = ^Ti(MtSt)at (5.60) 
The first equation is a Riccati equation and can be solved with a Cayley transform: 

Mt = {t-Mt){t + Mt)- 1 
which gives the linear equation 

Mt = StMt 

so we have, recalling that S t = F t Ff , 

M t = F t N 

Coming back to A4, we get 

M t = (t + Mo - F t (l - Mo)) (1 + Mo + F t (t - Mo))' 1 (5.61) 
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Let us now compute the phase term. We introduce 

x f = l + Mo±F t (l-M ) (5.62) 

Using the following properties 

XT = xt-2Ft(t-M ) (5.63) 

TiS t = (5.64) 

S t = F t F t - 1 (5.65) 



we have 
so we get 



Tr(M t S t ) = -2Tr {xtixt)' 1 )) (5-66) 



a 4 = a exp ^-^[logdetx+]o^ , a = l. 



(5.67) 



In formula (|5.67|) the log is defined by continuity, because we shall see that xt * s 
always non singular. 

Until now we just compute at the formal level. To make the argument rigorous 
we state some lemmas. 

It is convenient here to introduce the following notations. 

sp + (2n, C) is the set of complex, 2n x 2n matrices M. such that M. = JM where M 
is symmetric (that means that M. is a complex Hamiltonian matrix) and such that 
QM is positive non degenerate. 

Sp + (2n, C) is the set of complex, symplectic, 2n x 2n matrices M such that the 
quadratic form 

z*-+$l(l- ArW) z ■ Jz 
is positive and non degenerate on C 2n . 

Lemma 5.2 If F is a real symplectic matrix and M. E sp + (2n, C), then 1 + M. + 
F(l — A4) is invertible. 

Proof 

Write M = JM. It is enough to prove that the adjoint 1 + F T + M*(t - F T ) is 
injective. But M* = -MJ. So if z G C 2 " is such that (1 + F T + M*(l - F T ))z = 
then we get 

((1 + F T )z - MJ(1 - F T )z ■ J(l - F T )z) = (5.68) 

But, using that F is symplectic, we have that (1 + F)J{1 — F T ) = FJ — JF T is 
symmetric so taking the imaginary part in (|5.68|) . we have 

3 (MJ(1 - F T )z ■ J(l - F T )z) = 0. 

Then using that 3M is non degenerate, we get successively (1 — F T )z = 0, (1 + 
F T )z = and z = 0. □ 
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Lemma 5.3 Assume that -1 is not an eigenvalue of N '. Then M G Sp + (2n, C) if 



then we can easily see that M is symplectic if and only if Ai is Hamiltonian. 
Now, using A/* = Af* ,T , we get 

11 - ArW = 2(1 - M)-\M -M)(l + MY 1 

If z = (1 + M)z' we have, using (J(l + A4~)) T = - J(l - M), 

(1 - JT 1 M)z ■ Jz = 2(M - M)/ • I 7 . 

So the conclusion of the lemma follows easily from the last equality □ 

The last lemma has the following useful consequence. 
Let us start with some A^o £ Sp + (2n, C) without the eigenvalue -1. It not difficult 
to see that M.t z — ~z if and only if Ai u = —u, where u = (xt)^ 1 ?- In particular 
for every time t, -1 is not an eigenvalue for A4 t . 

Furthermore, using that Mt = F t Mo, we have J\f t Af t = A7 Ao- So we get that the 
matrix M. t G sp + (2n, C) at every time t. 

If M.q has the eigenvalue —1 it is no more possible to use the Cayley transform but 
we see that M.t is still defined by equation ([5.61)1 (from lemma (|5.2|) ) and solves the 
Riccati equation ()5.59|) . 

Now we want to discuss in more details the phase factor included in the term 
a t and to consider the limiting case M. Q = to compute the Weyl symbol of the 
propagator U®. So doing we shall recover the Mehlig- Wilkinson formula, including 
the phase correction Maslov-Conley-Zehnder index). 
Let us denote 



Lemma 5.4 Let us consider t i— > F t a path in Sp(n, R). Then for every Aio G 
sp + (n, C) we have, for the real part: 




Assuming that M + 1 is invertible, and 

M = (l-Af)(l + Af)-\ 




Hence we have 





18 



If F t is t '-periodic, then 

nT 5(F t ,M ) 



ds = 2inv 



l 5(F t ,M ) 

with v G Z. Furthermore, v is independent on M.q G sp + (2n, C) and depends only 
on the homotopy class of the closed path t i— > F t in Sp(2n). 

Proof 

For simplicity, let us denote 8(s) = 5{F S ,A4 ) and 

h{t) = [ 5(s)5(s)- 1 ds, g{t) = e- h{t) 5(t) 



g is clearly constant in time and g(0) = 5(0) = 1. Then we get $l(h(t)) = log \S(t)\. 
In the periodic case we have e h ^ = 1 so we have 

1 r5(F t ,M ) 

-ds = u, v G £. 



2nJ Q 5{F U M, 

By a simple continuity argument, we see that v is invariant by continuous deforma- 
tion on A4q and F t . □ 

We can now compute the Weyl symbol of R(F) when det(l + F) ^ 0. Let us 
consider first the case det(l + F) > 0. The case det(l + F) < is a little bit more 
complicated because the identity 11 is not in this component. 

We start with an arbitrary C 1 , path t i— > F t going from 1 (t — 0) to F (t — 1). 
It is known that Sp + (2n) = {F G Sp(2n), such that det(l + F) > 0} is an open 
connected subset of Sp(2n). So, we can choose a piecewise C 1 path F[ in Sp + (2n) 
going from 1 to F and A^o = is J- We have, using Lemma ()5.4|) . 



But det(l + F/) is never on [0, 1] and is real; so if e > is going to zero, the last 
term in r.h.s goes to and we get 

hmSl ffeM^I^,, 

£ \° \Jo S(F t ,ieJ) 

So we have proved for the Weyl symbol R(F, X) of R(F) the following Melhig- 
Wilkinson formula 

R(F, X) = e i7rI/ |det(ll + F)\~ 1/2 exp (-iJ(l - F)(l + F)~ 1 X ■ X) (5.70) 
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Let us now consider F G Sp_(2n) where Sp_(2n) = {F G Sp(2n), such that det(l+ 
F) < 0}. Here we shall replace the identity matrix by 

F ' ~ { "I 

for n = 1 and F° = F 2 ° <g> \2n-2 for n > 2 where \% n -% is the identity in IR 2n ~ 2 . 
Let us consider a path connecting 11 to F° then F° to iq = F. Because Sp(2n) is 
open and connected we can find a path in Sp(2n) going from F° to F± = F and that 
part does not contribute to the phase by the same argument as above. 
Let us consider the model case n = 1. The following formula gives an explicit path 
in Sp(2). 

p, _ ( costn — sintn \ / T](t) 

* y Smt7T COSt7T J y 

where rj(t) is analytic on a complex neighborhood of [0,1], 77(0) = 1, 77(1) = 2, 
1 < v(t) — 2 for f G [0,1]. A simple example is i](t) = 1 + 1. Then we can compute: 



Lemma 5.5 

""0 " \Jo S(F t ,ieJ) 



Um a I I 1 i[ F y £j } dt I = ^ (5.71 ! 



Proof 

Let us denote F' t = R(t)B(t), where R(t) is the rotation matrix of angle tir and for 
£G[0,1[, 

f £ {t) = det[(l - ieJ) + - ieJ)] 

We have 

fo(t) = det(R(-t) + B(t)) = 2 + cos(tTr) (l+t+ 1 



l + t 



fo has exactly one simple zero t\ on [0, 1], /o(£i) = 0, /o(ti) 7^ 0). 

This is easy to see by solving the equation cosi7r = h(t), for a suitable h, with a 

geometric argument . 

Then by a standard complex analysis argument (contour deformation) we get the 
equality (lo~7T| . □ 

So in this case we have the formula ()5.70|) for the contravariant Weyl symbol of 
R(F), with index v G Z + 1/2. Summing up the discussion of this section we have 
proved: 

Theorem 5.6 We can realize the metaplectic representation F t— > -R(-F) of the sym- 
plectic group Sp(2n,M) into t/ie unitary group of L 2 (M. n ) by taking for every F a C 1 - 
path 7 going from l(t = 0) to F(t = 1) and solving explicitly the corresponding 
quadratic Schrodinger equation for the Hamiltonian H t (z) = — l/2JF t Ff 1 z ■ z. So 
let us define R 1 {F t ), the propagator at time 1 obtained this way. 
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If ^ is another path going from t (t = 0) to F (t = 1) then there exists an index 
iV(7, 7') G Z such that 

Ry(F) = e i7rAf(7 ' 7,) ^(F) 

The metaplectic operator R(F) is the two valued unitary operator ±R 1 (F). 
Moreover if det(t + F) ^ 0, R(F) has a smooth contravariant Weyl symbol R(F, X), 
given by formula \5. 70Q , where v G Z if det(l + F) > and v 6 Z + 1/2 z/ 
det(l + F) < 0. 

It will be useful to translate the above theorem for the covariant Weyl symbol. 
Before that we start to discuss the general case, including det[F± 1] = 0. The Weyl 
symbol of R(F) may be singular, so it is easier to analyse it using coherent states 
(for our application it is exactly what we need). 
Let us recall that 

R*(F,X) = (2irh)- n [ (<p g+x \R(F)(p t )e-^< x ^dz (5.72) 

JR 2n 

Let us denote U\ = R(F) and U\ the contravariant Weyl symbol at time 1 con- 
structed above, such that Uq(X) = e _e ' x ' . 

For every e>0we have computed the following formula for the contravariant Weyl 
symbol: 

U £ {X) = a £ exp (iM £ X ■ X) , where 

M £ (X) = =-J(t + teJ-F(t-tsJ))(t + isJ + F(l-ieJ))- 1 (5.73) 
a £ = det(l + is J + F(l- is J))- 1 ' 2 (5.74) 

At the end we get the result by taking the limit: 

lim (<p g+x \Ul<Pz) = (<p z+ x\R(F)<p z ) (5.75) 

The computation uses the following formula for the Wigner function, W z , z+ x{Y), of 
the pair (y? 2 ,y5 2+ x)- 





Y-z-X 


- 






2 



ia{X,Y-t) 



W ZtZ+x {Y) = 2 2n ex V 
So, we have to compute the following Fourier-Gauss integral 
(ip z+x \U&z) = 2 2n (2n)- n a £ [ dY exp I iM £ Y ■ Y 
So we get 

{•P.+xPt'P.) = 2" (det(l - iM'))- 1 ' 2 a'. 



Y-.-X 
2 



(5.76) 



ia(X,Y-±) 
(5.77) 

(5.78) 
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x exp 



z + 



X 



1 



ia{X, z) + (l-iM e 



X-iJX 



z + 



X-iJX^ 



2 V ' /V /v 2 /v 2 

Now we can compute the limit when e \ 0. We have 

11 - iM £ = 11 + iJ(l + zeJ - F(l - ieJ)){t + ieJ - F(l - ieJ))- 1 (5.79) 
Then we get, using that (1 + F + i J(t — F)) is invertible (see Lemma f5.2j) . 

lim(ll - iM £ )~ l = (1 + F)(1 + F + U{t - F))' 1 (5.80) 

and 



limdet(l-iM £ )- 1/ V = (det(l + F + i J(l - F)) 



-1/2 



(5.81) 



So, finally, we have proved the following 



Proposition 5.7 The matrix elements of R(F) on coherent states (p z , are given by 
the following formula: 



(<p z+x \R(F)<p z ) = 2 n (det(l + F + U{1 - F)y 



-1/2 



x exp 
where 





X 


- 


2 H 




2 



1 . . X -iJX s . X -iJX s . 
+ + )-{z + )) (5.82) 



K F = (1 + F)(l + F + zJ(l - F))" 1 



(5.83) 



Now we can compute the distribution covariant symbol of R(F) by plugging formula 

in formula (l5~72l . 
Let us begin with the regular case det(l — F) ^ 0. 

Corollary 5.8 7/det(l — F) ^ 0, the covariant Weyl symbol of R^{F) is computed 
by the formula: 

R*(F,z) =e mM |det(ll-F)|- 1/2 exp (~ J(l + F)(l - F)~ 1 z ■ z\ (5.84) 

where fi = v+ ^,b , &Zisan index computed below in formulas \5.88\) . i5.89\) . 
Proof 

Using Proposition (J5.7|) and formula (|5.72|) . we have to compute a Gaussian integral 
with a complex, quadratic, non degenerate covariance matrix (see 
This covariance matrix is Kp — \ and we have clearly 



Kp — 1 = —iJ(t — F)(t + F + i.J{l — F))^ 1 = -(1-iA)" 1 
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where A = (1 + F)(t — F) 1 J is a real symmetric matrix. So we have 

&(K F -1) = -(fl + A 2 )" 1 , $s(K F - 1) = -A(l + A 2 ) -1 (5.85) 

So that t — Kp is in the Siegel space S 2n and Theorem (7.6.1) of [TT] can be applied. 
The only serious problem is to compute the index //. 

Let us define a path of 2n x 2n symplectic matrices as follows: Gt = e t7Tj2n if 
det(l - F) > and: 

G t = G\ g> e^ j2 "~ 2 if det(l - F) < 0, where 



G 



( v(t) 
V W) 



where rj is a smooth function on [0, 1] such that 77(0) = 1, rj(t) > 1 on ]0, 1] and 
where Jm is the 2n x 2n matrix defining the symplectic matrix on the Euclidean 
space R 2n . 

G\ and F are in the same connected component of Sp^.(2n) where Sp Vr (2n) = {F G 
Sp(2n), det(l — F) 7^ 0}. So we can consider a path s 1— > F s ' in Sp 7lr (2n) such that 
Fq = Gi and F[ = F. 

Let us consider the following "argument of determinant" functions for families of 
complex matrices. 

0[F t ] = arg c [det(l + F t + U{1 - F t )) (5.86) 
p[F] = arg + [det(l- Kp)- 1 } (5.87) 

where arg c means that t 1— > 0[F t ] is continuous in t and 8[t] = (Fq = t), and 
S t— > arg_ | _[det(5')] is the analytic determination defined on the Siegel space S 2n such 
that arg + [det(S')] = if S is real (see [TT], vol.1, section (3.4)). 
With these notations we have 

,.m=m. (5.88) 

Let us consider first the case det(l — F) > 0. 

Using that J has the spectrum ±i, we get: det(l + G t + iJ(l — G t )) = A n e ntm and 
t-K Gl = t. 

Let us remark that det(l - K F )~ l = det(l - F)~ 1 det(ll - F + iJ(l + F)). Let us 
introduce A(E,M) = det(l-E + M(l + E)) for E G Sp(2n) and .M G sp + (2n,C). 
Let consider the closed path C in Sp(2n) defined by adding {G t }o<t<i and {F^} < s <i. 
We denote by 2iru the variation of the argument for A(«, Ai) along C. Then we get 
easily 

p(F) = 9[F] +2vrz/ + mr, n G Z. (5.89) 

When det(l — F) < 0, by an explicit computation, we find arg + [det(l — F Gl )] = 0. 
So we can conclude as above. □ 
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Assume now that F has the eigenvalue 1 with some multiplicity 2d. We want to 
compute the temperate distribution B}(F) as a limit of B}(F e ) where det(l — F e ) ^ 
0, \/e > 0. 

Let us introduce the generalized eigenspace £' = (Jker(l - F) j (dim[£']=2d). and 

£" its symplectic orthogonal in M? n . We denote F' the restriction of F to £' and 
F" the restriction to £". We also denote by J' and J" the symplectic applications 
defined by the restrictions of the symplectic form a: a(u, v) = J'u ■ v, Viz, v G £' and 
the same for J". Let us introduce the Hamiltonian maps V = (1 — + 
and L'' £ = V — eJ'. 

It is clear that det(Z/' e — 1) ^ for < £ small enough, so we can define 
F'' £ = (1 + L'' e )(l - L /,e ) _1 . Let us remark that 

Q'- £ := J'(F'' e - l)-^^' 6 + 1) = -(L'J' + e)- 1 (5.90) 

is a symmetric non degenerate matrix in £' defined for every e > 0. 

Lemma 5.9 We have the following properties. 

1) F' ,£ is symplectic. 

2) limF'- 6 = F'. 

3) Fore^O, small enough, det(F'' £ - 1) ^ 0. 
Proof. 

1) comes from the fact that L' ,e is Hamiltonian. 

2) is clear. 

For 3), let us assume that F'' £ u = u. Then we have V' £ u = hence J'L'u = eu. 

Now, choose < e < dist{0, spec( J'L')\{0}} then we have u = 0. 

Finally we define F £ = F l,£ ® F" ,£ . It is clear that F £ satisfies also properties 1), 2), 

3) of the above lemma with F £ in place of F' )£ . □ 

Proposition 5.10 Under the above assumption, the covariant symbol of R(F) has 
the following form: 

R*(F, Zl , z 2 , z") = e"^|det(Il - F")\- l l 2 5( Zl ) 
x exp + F)(F- t)-\z 2 + z") ■ (z 2 + z")\ (5.91) 

where z := ((zi, z 2 ), z") is the decomposition of the phase-space for which F = 
F' ® F" and z' — z\ + z 2 , with z 2 e lm(F' — 1), z\ G Im(F' — l)- 1 , the orthogonal 
complement in £' for the Euclidean scalar product. 5{z\) denotes the Dirac mass at 
point z\ = 0. /ii G Z + 1/2 is given as follows: ll\ = fx" + sg .^ where fx" is the limit 
of the fx index for F £ , computed in \5.8J$ , and sg + Q' is the limit for e — > of the 
signature of Q'' £ , defined in A5.9U\) . 



24 



Proof 

We use the same kind of computation as for Corollary (|5.8j) . The new factor comes 
from the contribution of £'. So we can forget £". So we assume that £' = M 2n and 
and we forget the superscripts '. We have to compute the limit in the distribution 
sense of F$(F £ ) which was computed in Corollary (|5.8|) . 

Let us consider a test function / £ iS(M 2n ) and its Fourier transform /. Using 
Plancherel formula, we have 

/ exp ( -7.7(1 + F £ )(t - F £ )- l z ■ z) f(z)dz = 
7r2™ V 4 / 

2 n (7r) 2n |det(g e )|- 1 / 2 e sg(Q£) ^ I exp + F £ )-\l - F £ )J( ■ () f(()d( (5.92) 

So we get the result by taking the limit for e — > in (|5.92j) . 

Let us remark that we have used that J'(ll + F)(t — F) _1 is an isomorphism from 
ker(l - F') 1 - onto itself. □ 

In our paper [2] the leading term for the return probability and for fidelity on 
coherent states is computed with a Gaussian exponential defined by the quadratic 
form which was defined in Proposition 15.71 

lF (X) := i (K F (1 - U)X ■ (1 - U)X - \X\ 2 ) (5.93) 

In our application we shall have X = z t — z where 1 1— > z t is the classical path starting 
from z. So we need to estimate the argument in the exponent of formula (|5.82j) for 
z = 0. 



Lemma 5.11 we have, VX 6 



s2n 



M( lF (X)) < - \ X } 2 (5.94) 
2(1 + s F ) 

where s F is the largest eigenvalue of FF T (F T is the transposed matrix of F) . 
Proof 

Let us begin by assuming that det(l + F) 7^ 0. Then we have 

K F = (1 + IN)' 1 , where N — J(l — F)(l + F)' 1 . 
So we can compute 

&{K F ) = {t + N 2 )- 1 = K F K F and Q(K F ) = -N(l + N 2 )' 1 . 
So, we get, 

lF (X) = - ((1 + JN)K F K F {1 - NJ)X ■ X - 2\X\ 2 ) (5.95) 
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By definition of Kp, we have 

(1 + JN)K F = 2 ((1 + i^F' 1 + 1-U) ^ (5.96) 
Le us denote Tp = ((1 + iJ)F _1 + 11 — iJ)^ 1 . We have, using that F is symplectic, 

T F '~ l Tp l = 2{F~ 1 ' T F~ 1 + 1). 

hence we get 

TpT F = {2{F~ 1 ' T F- 1 + I))" 1 

and the conclusion of the lemma follows for det(l— F) ^ hence for every symplectic 
matrix F by continuity. □ 
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